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Abstract 

We discuss the Drell-Yan process at a measured transverse-momentum Qt of the 
produced lepton pair in collisions of transversely polarized protons and antiprotons, 
to be observed at the proposed spin experiments at GSI. The large logarithmic con- 
tributions from multiple soft gluon emission, accompanying the Drell-Yan mechanism 
at small Qt, are resummed to all orders in QCD perturbation theory up to next- 
to-leading logarithmic (NLL) accuracy. Numerical evaluation shows the impact of 
the NLL as well as LL effect on the dilepton Qt spectra. For the corresponding 
Qr-dependent spin asymmetry £^tt{Qt)-, the LL effect gives significant modification 
while the NLL effect is marginal, leading to QCD prediction that £^tt{Qt) at GSI is 
flat at small and moderate Qt and almost equals the conventional asymmetry Att 
associated with the QT-integrated cross sections. This flat behavior in turn allows 
us to use analytic saddle-point evaluation of the resummation formula in the limit 
Qt 0, not only to obtain quantitative estimate of ^tt{Qt), but also to clarify 
mechanisms behind the relation ^tt{Qt) ~ Att characteristic of pp collisions at 
GSL 
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Recently there has been much progress to unveil the last unknown parton distribution 
of nucleon at the leading twist, the transversity distribution 5q{x) [Il[2]: for example, the 
calculation of the lowest two moments of 5q{x) by lattice QCD simulation is updated [3], 
and the first global fit of 5q{x) is given [1] using the semi-inclusive deep inelastic scattering 
(SIDIS) data, in combination with the e^e~ data for the associated (Collins) fragmentation 
function. These results indicate that the transversity distributions for u and d quarks are 
sizable, but are fairly small compared with their positivity bound (Soffer bound) [5]; still, 
these results are subject to certain assumptions and uncertainties, and it is important to 
proceed further toward the determination of the transversity. 

As is well-known [HIS], the transversely polarized Drell-Yan (tDY) process provides an- 
other promising way to access 5q{x). The measurement of the tDY cross section is proposed 
in the future experiments at GSI [6], where the pp collisions at moderate energy mainly 
probes the products of two quark transversity-distributions, 5q{xi)5q{x2)i in the "valence 
region". The corresponding double-spin asymmetries estimated at the leading order (LO) 
in QCD are large enough to be measured at GSI [7], in contrast to the complementary case 
of tDY in pp collisions at RHIC, where the asymmetries are predicted to be rather small be- 
cause 5q{xi)5q{x2) + (1^2) is probed in the "sea-quark region" (see [8l[9l[10]). As for the 
higher-order QCD corrections for tDY in pp collisions, the NLO corrections have been stud- 
ied recently [11]; also, the "threshold resummation" has been applied in order to sum up, 
to all orders in a^, the soft-gluon emission contributions that are logarithmically enhanced 
near the threshold of the partonic process [12]. The effects of these QCD corrections are 
small for the double transverse-spin asymmetries in pp collisions at the kinematical regions 
corresponding to the GSI experiments [HI [12], suggesting that the large LO asymmetries 
obtained in [7] are rather robust. 

All the above previous studies of tDY in pp collisions considered the case in which the 
transverse-momentum Qt of the produced lepton pair are unobserved. Experimentally, 
however, the bulk of events is produced in the small Qt region. Therefore, it is desirable 
to develop theoretical predictions of tDY at a measured Qt in pp collisions for the de- 
tailed comparison with the data in the future GSI experiments. When Qt <^ Q with Q 
the dilepton mass, the fixed-order perturbation theory breaks down due to the appearance 
of large logarithms hi{Q'^ /Q\) multiplying a^, and we have to deal with the relevant DY 
cross sections at small Qt in an all-order resummation framework in QCD. The corre- 
sponding "Qr-resummation" has formally some resemblance to the threshold resummation 
mentioned above, but represents other contributions associated with different "edge region" 
of phase space. The Q^-resummation for tDY has been formulated recently by the present 
authors [131 El [ID], summing the corresponding large logarithms up to next-to- leading log- 
arithmic (NLL) accuracy, in the context of the pp collisions at RHIC and J-PARC. There 
we demonstrated that profound modifications arise in the DY production at small Qt, 
driven by partonic mechanism that induces the large logarithmic contributions; in par- 
ticular, the interplay between the Sudakov factor resumming multiple soft gluon emission 
and the DGLAP evolutions of parton distributions yields the "amplification" of the double 
transverse-spin asymmetries in the small Qt region at RHIC as well as J-PARC kinemat- 
ics, resulting in their values larger than the conventional (fixed-order) NLO asjTiimetries 
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[8] calculated in the Qr-unobserved case. In this Letter we apply our Qr-resummation 
formalism to the case of pp collisions, and study the corresponding double transverse-spin 
asymmetries in the small Qt region for tDY foreseen at GSI. 

Detailed derivation of the tDY cross sections in our Qt resummation formalism is given 
in [in] for collisions of spin-1/2 hadrons hi and h2, i.e., hih2 — > l~^l~X, and the correspond- 
ing results for pp collisions can be obtained by trivial substitutions. The spin-dependent 
{ATda = {da^^ — da^^)/2) and spin-independent {da = {da^^ + da^^)/2) parts of the dif- 
ferential cross section are expressed as 



(At) da /oj,\ /o\ 



2 



cos(20)/2) 



dQ'^dQldyd<p ' '?,N,SQ'^ 



{At) X'^^\Q't. Q\ y) + (At) ^ (Q^, Q\ y) 



Here and below we follow the convention of [21 |TU] for the basic quantities entering the 
formulae: \fS and y are the total energy and dilepton's rapidity in the proton-antiproton 
CM system, and the prefactor cos(2(/))/2 specific to the spin-dependent part shows the 
characteristic dependence [1] on the azimuthal angle of one of the outgoing leptons with 
respect to the incoming nucleon's spin axis. The first term {AtX^^^ or X^^^) gives the 
dominant contribution when Qt ^ Q-, containing all the logarithmically-enhanced contri- 
butions a^\n^{Q'^ / Q"]^) / Qt and, at the NLL accuracy, has to be evaluated by resumming 
the first three towers (m = 2n — 1, 2n — 2, 2n — 3) of these large logarithmic contributions 
to all orders in as- The second term {AtY or Y) is free of such contributions, and can be 
computed by fixed-order truncation of the perturbation theory. 

The NLL resummed component {At)X^^^ is obtained through various kinds of elabo- 
ration [131 El [ini [HI [IS] of the Collins-Soper-Sterman (CSS) resummation formalism 
Introducing the impact parameter h space, which is conjugate to the Qt space, we hav^ 



as{Q^) 



(3) 



where the DY scaling variables are denoted as = 

JoibQr) is a Bessel function, and 60 = '2e~'^'^ with the Euler constant. 6H, H, and K 
denote the products of the NLO parton distributions of proton and antiproton, summed 
over the massless quark flavors q with their charge squared e^, as 



In this paper, we set i^r — hf = Q for the rcnomiahzation and factorization scales /ir and /ip- 
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and K{xi,X2;n^) = Y.q^l fi'(a;i, /i^) [g(x2, /i^) + g(x2, /i^)], where 6q{x,fi^), q{x,fi^), and 
gi^x^jji^) are the quark transversity, quark density, and gluon density distributions of the 
proton, respectively; note that there is no transversely-polarized gluon distribution at the 
leading twist. The Sudakov factor e^'^^'^^ and the coefficient functions {AT)C^j\z) are 
perturbatively calculable, and the former gives the all-order resummation of logarithmically 
enhanced contributions due to multiple emission of soft and/or coUinear gluons from the 
incoming partons. The exponent S{b, Q) is expressed as 

S{b,Q) = ^—h^'\X) + hW{X), (5) 

where the first and second terms collect the LL and NLL contributions, respectively, in 
terms of the two functions h^'^\\) and h^^\\), which depend on as{Q^) only through 

A = PoasiQ') HQ'bVbl + 1)= Poas{Q')L , (6) 

with (3q the first coefficient of the QCD f3 function. When Qx <^ Q, L plays the role of the 
large-logarithmic expansion parameter in the b space, as 6 ~ ^/Qt- In this relevant region, 
A can be as large as 1 even for ^^(Q^) <^ 1, and the ratio of two terms in ([5]) is of 0{as)] note 
that the NNLL or higher-level corrections, which are down by as or more, are neglected in 
([5]) [H [To]. The Sudakov exponent is independent of process as well as scheme to the NLL 
accuracy [IH [HI [17] , so that ([5]) gives the universal Sudakov factor common to ([2]) and ([3]). 
The explicit form of h^^^X) and h^^\X), as well as the coefficient functions {AT)Clj\z) 
in the MS scheme, can be found in [131 El EH]- It is implicit in ([2]) and ([3]) that the b 
dependence of SH{xi,X2]bl/b'^), H{xi, X2;bl/b'^) and K{xi,X2]bl/b'^), associated with the 
NLO perturbative evolution of the parton distributions from the factorization scale fip = Q 
to the scale bo/b, is also organized in terms of ([6]) to ensure the consistent NLL accuracy; 
i.e., the customary NLO evolution operators for the distributions are expanded up to the 
NLL term, where the LL term proves to be absent P, [ID]. The Fourier transformation to 
the Qt space in ([2]), ([3D is performed along a contour C in the complex b space P [TU l [T ^ [T^ . 
avoiding the singularity of the Sudakov exponent ([S]) at A = 1, which is associated with 
the Landau pole in the perturbative running coupling. This singularity in the b space 
signals the onset of additional nonperturbative phenomena at very large values of \b\, and 
the corresponding nonperturbative effects are complemented in ([2D, ([3D by introducing 
a Gaussian smearing function exp{—gf^pb'^) with a parameter gNP, following the usual 
procedure [131 [HI [13 [IS]- This may be interpreted as representing "intrinsic transverse 
momentum" of partons inside proton, and we use the same smearing function for both 
polarized and unpolarized cases, following our previous works [9l [10] for tDY in pp collisions. 

The "regular component" {At)Y is determined by the matching procedure expanding 
(121), ([3D in powers of as(Q^) and assuming gj^p — s> in perturbation theory, so that ([ID 
coincides exactly with the fixed-order result of the corresponding polarized and unpolar- 
ized differential cross sections, up to 0{as) [I3l [T5l [H] . Namely, the LO cross section for 
Qt > 0, {AT)da^'-' / dQ'^dQ^dyd(f), which is of 0{as) because the finite Qt of the lepton 
pair is provided by the recoil from the gluon radiation, is given by ([ID with the replace- 
ment [At)X^^^ {At)X^^^\-fo, where {At)X^^^\-po denotes the terms resulting from 
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the expansion of the resummed expression up to the fixed-order ^^(Q^). According to 
the structure of ([1]) via the matching with the LO cross sections, we refer to ([1]) as the 
"NLL+LO" prediction, and this gives the tDY differential cross sections in the MS scheme, 
which are well-defined over the entire range of Qt- 

The ratio from ([T]) yields the double transverse-spin asymmetry in tDY as 

2 X''^^(Ql,Q\v) + Y(Ql,Q\y) 

for measured Qt, Qi y, and (p. To the fixed-order as without the soft gluon resummation, 
(At) X^^^ (Aj-) X^'"'"|fo as discussed above, and ([7]) reduces to the LO prediction 
■^tt{Qt) for Qt > 0. We also introduce the asymmetry in terms of the NLL resummed 
components, ([2]) and ([3D, which, respectively, dominate the numerator and denominator in 
([ZD when Qt < Q: 



yNLL 



W.) = 1 00.(2^) ^^^""<f''^''^'. (8) 



As emphasized in [9l [TO] in the context of the pp collisions, the Qt limit of ([2D and ([3D 
deserves special attention: at Qt = 0, the b integral of ([2D and ([3D is controlled by a saddle 
point and can be evaluated analytically as 



(AT)xNLL(o,g^l/) 



/ 2^ ^-C^'>H\sp)+hW(Xsp) 



^Q^f3oas{Q^)]l ((oy'iXsp) 

(9) 

where C^°\X) = -\/[poas{Q^)] - h^^\\) / as{Q^) + [gNphl/Q'^]e^'^^^^^^^ '^\ and the saddle- 
point value of ([SDQ, Xsp = Poas{Q^) ln{Q%lp/bl), is defined by the condition C^°^'(A5p) = 0. 
The saddle-point formula ([ID is exact up to the 0{as) corrections that actually correspond to 
the NNLL contributions in the region Qt ~ (see the discussion below (Q, and also [9l [TO] 
for the details). Note that the gluon distribution completely decouples from X^^^{0, Q"^, y) 
at the NLL accuracy for Qt ~ 0. The prefactor inside the square bracket in the RHS 
involves "large perturbative effects" due to the Sudakov factor, as well as the Gaussian 
smearing factor with gNP'-, the former contribution drives the well-known asymptotic be- 
havior [18] of the DY cross sections, ~ (A^cd/Q^)"^"^^^^^"^ with a = 4 /(2vr/?o), for 
Q ^ Aqcd- Because this prefactor is common to both the polarized and unpolarized cross 
sections, we obtain the remarkably compact formula for the Qt limit of ([HD [HI fTO] : 

(Q^ _ 0) - - coK20)-^^^p^-^p^ . (10) 

Using the formulae described above, we study the behavior of tDY to be observed in pp 
collisions at GSI. The PAX Collaboration has proposed the tDY experiments in pp collisions 



t For the kinematics of our interest, we have the saddle point well above 6 = 0, and we can use the 
definition A = f3oas{Q^)ln{Q^b'^ /b^), up to the exponentially suppressed corrections to ((Hj) (see [TU]). 
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Figure 1: The transversity distributions for (a) w-quark and (b) d-quark. The sohd and 
dashed curves plot the NLO distributions corresponding to the Soffer bound and the relation 
(fTTj) . respectively. The dot-dashed curve shows the LO distribution corresponding to the 
Soffer bound and the shaded area shows the result of the LO global fit in . 

at S" = 30 and 45 GeV^ in the fixed-target mode, and those up to 5 = 210 GeV^ in the 
collider mode [B]. Those GSI-PAX experiments will probe 0.2 < Q/\/S < 0.7, and thus 
the transversities in the "valence region" in a wide range of x. To compute the tDY cross 
sections ([1]) at the NLL-I-LO accuracy with these GSI kinematics, we have to specify the 
NLO parton distributions to be substituted. We use the NLO GRV98 distributions [19] 
for the unpolarized quark and gluon distributions /i^) and g{x,^'^). For the NLO 
transversity distributions 5q{x, iJ?), we consider the two typical assumptions that have been 
used in the literature [71 El [HI El [IDl [IH : at a low input scale /xo < 1 GeV, these 
assume the saturation of Soffer's inequality [5] as Sq{x,fil) = [q{x,fil) + Aq{x, ^l)]/2, and 
the relation, 

(5g(x,/io) = Ag(x,/io) , (H) 

exact in the non-relativistic limit, respectively, and their QCD evolution from fiQ to a higher 
scale /i is controlled by the NLO DGLAP kernel [2T] for the transversity; here Aq{x, fi"^) 
denote the longitudinally polarized quark distributions. The first case yields 5q{x,ii^) that 
satisfies Soffer's inequalitjH and provides an upper bound on the transversities [8]. The 
second case (fTTj) is suggested also by the estimates from relativistic quark models for nu- 
cleon [21 [71 [22|. For the input functions in the RHS of these two assumptions, we take 
the NLO GRV98 distributions g(x, /Iq), as noted above, and GRSV2000 ("standard sce- 
nario") distributions Ag(x, /^q) [20] with /Xq = 0.40 GeV^. The obtained NLO transversity 
distributions for u and d quarks, and x6d{x, fi"^), are shown in Figs. 1 (a) and 

(b), respectively, as a function of x with /i^ = 2.4 GeV^. The solid curve shows the re- 
sult corresponding to the Soffer bound, and the dashed curve shows the result using f[TT]) . 
For comparison, we also depict the shaded area which represents the result (a one-sigma 

^ The obtained (5g(x,/L<^) is actually very close to + Aq{x, fi^)]/2, except for small x (< 0.2). 



6 



confidence interval) for tlie LO transversity distributions extracted tfirougfi tlie global fit 
to tlie data j^, and the dot-dashed curve which shows the LO transversity distributions 
corresponding to the Soffer bound with the LO inputs of [191 ISO]- We note that the above 
assumption for the Soffer bound, Sq{x,fil) = [q{x,fil) + Aq{x, fil)]/2, following the litera- 
ture 0, [131 El [ini [m [12] yields the positive polarization for the d quark, 6d{x, iJ?) > 0, in 
contrast to the result using ([TTI) . but the sign of the polarization cannot be detected in tDY 
for pp collision (see (jl])). For convenience of comparison, the solid and dot-dashed curves 
in Fig. 1(b) show the results multiplied by —1. 

The results in Figs. 1 (a) and (b) indicate that the empirical LO transversities are 
smaller compared with the transversities corresponding to the Soffer bound, in particular 
for the u quark. On the other hand, the NLO transversities using f[TT]) lie slightly outside 
the one-sigma error bounds of the empirical fit. The two NLO sets as well as the empirical 
fit have {Su{x, yU^))^ ^ {5d{x, /i^))^ in the valence region relevant at GSI kinematics, so that 

SH{xi, X2\ /i^) ^ e^(5u(xi, ft^)5u{x2, fj.^) , (12) 

for (jlj). Hence the cross section ([T|) and the asymmetry ([7]) at GSI allow a direct access to 
\6u{x, fi'^)\. The NLO transversities using f[TTl) satisfy Soffer's inequality for the u quark, 
but violate it for the d quark by a small amount, see the last footnote and Figs. 1 (a), (b)|l. 
However, this violation of Soffer's inequality will be harmless to our numerical estimates of 
the cross sections and asymmetries because of the dominance of the w-quark distribution 
noted above. We also mention that, strictly speaking, Soffer's inequality for the NLO 
distributions receives the additional scheme- dependent radiative corrections [5], and the 
corresponding corrections would modify the solid curve in Figs. 1 (a), (b) by a certain 
amount of 0{as)- 

In all the following numerical evaluation, we choose = for the azimuthal angle of a 
lepton; extension to other (p is straightforward by taking into account the cos(2</)) depen- 
dence displayed in the relevant formulae ([T|), ([7]), etc. For the nonperturbative parameter 
in (I2D and di, we use qnp = 0.5 GeV^ as used for the RHIC and J-PARC cases [131 [HI [ID]. 

Probing the transversity distributions through the tDY in the collider mode is more 
promising than that in the fixed-target mode, because at higher energies the description 
based on QCD perturbation theory is supposed to be more accurate [12]. In Fig. 2(a), we 
show the Qt spectrum of the dilepton in tDY with S* = 210 GeV^, Q = 4 GeV, and y = 0, 
corresponding to a kinematics in the collider mode. We use the transversity distributions 
corresponding to the Soffer bound (the solid curve in Fig. 1), which were used in similar 
study of the Qt spectrum at RHIC [131 [3 [IH]- The solid line shows the spin-dependent 
part of the NLL-I-LO differential cross section ([T|), multiplied by 2Qt, and the dot-dashed 
line plots the contribution from the NLL resummed part AtX^^^ given in ([2|). The dotted 
line plots the contribution from the LL resummed part. 



ATXLL(Q|,Q^^/) 



6H{xlx'i;Q') , (13) 



The corresponding u-, d-, and s-quark distributions satisfy Soffer's inequality. 
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Figure 2: The tDY with GSI kinematics, S = 2W GeV^ Q = 4 GeV, ?/ = and = 0, 
and with g^p = 0.5 GeV^, using the NLO transversity distributions which correspond to 
the Soffer bound: (a) The spin- dependent part of the cross section, Axda /dQ'^dQTdydcp. 
(b) The double transverse-spin asymmetries obtained using each curve in (a). 



which is obtained from by omitting the NLL terms with the nonperturbative inputs {qnp 
and transversity distributions) kept intact; note that the b dependence of the parton distri- 
butions in ([2]) is associated with the NLL-level terms, as mentioned below ([6]) [9l[T0l[T3]. The 
dashed line shows the Qt spectrum using the LO cross section Aj^da^^^ /dQ'^dQ'^dyd(f) intro- 
duced above which is divergent as Qt —>■ due to the singular terms oc \n{Q'^/Q^)/Q'^, 
oc 1/Qt- By resumming the singular large logarithms to all orders in Og, the Qt spectra are 
completely redistributed and well behaved, forming a peak at Qt ~ 1 GeV. In fact, around 
the peak region, the NLL+LO cross section is dominated by the contribution from the NLL 
resummed part AyX^^^. It is also remarkable that the NLL result is considerably enhanced 
compared with the LL result, though the integrations of these two results over Qt coincide, 
using A = at 6 = (see and P,[in]), up to the O (^^(Q^)) corrections associated with 
the coefficient function ATCql\ We note that the pattern similar to Fig. 2(a) is observed 
also in the corresponding Qt spectra for the unpolarized differential cross sections. 

The behavior of the NLL-I-LO cross section at Qt — 3 GeV in Fig. 2(a) is still affected 
by the "broadening" due to the nonperturbative smearing e~^^^^^ in ([2]); as a result, the 
NLL+LO cross section is larger than the LO one at Qt — 3 GeV, although the former was 
matched to the latter in the region where the logarithm ln(Q^/Q|,) is not large (see the 
discussion above ([7])). Because the separation between the "peak region" Qt ~ 1 GeV and 
the "matching region" Qt ~ Q is not so large for the present kinematics, the value Qt — 
3 GeV actually corresponds to the boundary between "smeared" and purely perturbative 
regimes. In fact, for higher and Q such that the peak and matching regions are largely 
separated, the NLL+LO cross section around Qt — Q is independent of the nonperturbative 
smearing and reduces to the LO one (see e.g. the results for RHIC in [131 El [ID])- Another 
consequence due to the moderate-energy kinematics is that the perturbative contributions of 
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and higher, involved in the NLL resummed component ([2]), are not neghgible around 
Qi- ~ Q at the quantitative level, as as{Q'^)/'n: ~ 0.07. It is worth noting that, even in much 
higher energy processes such as Z boson [15] and Higgs boson productions [T^ , the higher- 
order perturbative contributions in the resummed component remain numerically sizable 
also at rather large Qt well above the peak region. In the present case, however, the impact 
of those higher-order contributions can be even more significant, because the LO cross 
section at Qt ~ Q is very small (see Fig. 2(a)); it is actually much smaller than its canonical 
size, [a^ /?>NcSQ'^] x [a^i^Q'^) /txQ] ~ lO^^ pb/GeV^ (see (II])), by the extra factor ~ 0.01; here 
the extra small factor comes from the fact that, for the moderate energy 75 = 14.5 GeV, the 
2 — >■ 3 partonic processes emitting a real gluon with the "high" transverse momentum kx ~ 
Q are possible only from the initial-state partons with the large momentum fractions Xi^2 as 
X1X2 ^ 0.3 (^ Q'^/S), for which the products of the parton distributions in (jl]) are strongly 
suppressed (see Fig. 1). These points suggest that the accurate quantitative description 
around Qt ~ Q at GSI would eventually require the matching procedure taking into account 
also the O(a^) or higher contributions in perturbation theory, which is beyond our NLL+LO 
framework. One may expect that the mechanisms observed for the large Qt region in the 
"NLL-I-NLO" cross section of Z boson production [15] and in the "NNLL-I-NLO" cross 
section of Higgs boson production [H] could also provide a better treatment for the matching 
region Qt ~ Q in the present case. Nevertheless, in this paper, we insist on the NLL+LO 
framework because the tDY Qr-differential cross section is not known at NLO, and we 
restrict our quantitative discussion below the matching region. 

Figure 2(b) shows the asymmetries as functions of Qt, obtained by taking the ratio of 
each curve in Fig. 2 (a) to the corresponding Qt spectra for the unpolarized cross sections; 
i.e., the solid, dot-dashed, and dashed curves plot the NLL+LO {^tt{Qt) of ([7])), NLL 
{s^tt^^Qt) of dH])), and LO {^tt^Qt) below (I7D) asymmetries, respectively, and the dotted 
curve shows the LL asymmetry, obtained by using ( |T3l) and the similar formula for X^^, as 

-^TT = -cos(20) — = -cos(20) — — „ , (14) 

2 ^ X^HQ'T,Q^,y) 2 ^ i7(x?, x^; Q^) ^ ' 

which is constant in Qt- All asymmetries in Fig. 2(b), except the LO result, show similar 
fiat behavior with almost the same value ~ 25%. The results observed in Fig. 2(a) imply 
that these fiat behaviors are governed by the soft gluon resummation contributions. In 
particular, we have, for Qt around the peak region in Fig. 2(a), 

^tt{Qt) ^ ^tt'-^Qt) ^ ■<t'^^(0) , (15) 

and the fiat behavior, represented by the second (approximate) equality, refiects the fact 
that the soft gluon effects resummed into the Sudakov factor e^^'^'^^ are universal to the NLL 
accuracy between AtX^^^ and X^^^ in ([H]). In fact, similar fiat behavior is observed in pp 
collisions at RHIC and J-PARC kinematics [SI [ID]; note, remarkably, the present result for 
pp collisions turns out to be even fiatter. As Qt — ^ 0, away from the peak region of the cross 
section, ^tt{Qt) decreases slightly due to the terms cx: ln(Q^/Qf.)/Q^ contained in AtY 
and y in d?]). The difference between ^tt{Qt) and the LO asymmetry at Qt — 3 GeV 
reflects the above-mentioned discrepancy between the corresponding cross sections. 
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Figure 3: The ratio of the transversity to unpolarized quark distribution for w-quark at 
different scales /i. The upper and lower curves are obtained using the NLO transversity 
distributions corresponding to the Soffer bound and the relation (fTTj) . respectively. 

We here recall that, in the pp-collision cases at RHIC and J- PARC using the same 
input transversity distributions as in Fig. 2(b), the significant enhancement of ^tt^{Qt) 
(2^ £^tt{Qt)) compared with s^^^ has been found [9l[T0]. Such enhancement is not seen in 
the present pp-coUision case. This fact indicates that the large NLL-level corrections shown 
in Fig. 2(a) are canceled in the asymmetry ■s^tt^{Qt) with the corresponding corrections 
to the unpolarized cross section, although this is not the case in the pp collisions. 

To clarify the reason behind this remarkable difference between the pp- and pp-coUision 
cases, the saddle-point formula (ITOl) is useful. Similarly to the pp-collision case P ITO] . 
the property f|T5|) allows us to use ffTOj) as a sufficiently accurate estimation of £^tt{Qt), 
£^tt^{Qt) (see Fig. 4 and Table 1 below). We obtain ho/hsp = 1.0 GeV for the scale of the 
transversity and unpolarized distributions in the numerator and denominator of ffTOj) at the 
kinematics of Fig. 2, using the condition that was noted below (Efl. £^tt of dH is different 
from flTU]) . only in the scale Q {= 4 GeV for Fig. 2) of the parton distributions. There- 
fore, the possible enhancement of s^tt{Qt)i ■^tt^{Qt) in comparison with s^^^ can be 
understood as a result of QCD evolution of the parton distributions from bo/bsp to Q, with 

^ b^/b'^gp: in fact, the corresponding enhancement in pp collisions at RHIC and J-PARC 
arises because of very different behavior of the sea-quark components under the evolution 
between transversity and unpolarized distributions ^ [TU] ; note that the transversity distri- 
butions obey a non-singlet-type evolution even for the sea-quark components because there 
is no gluon transversity distribution. On the other hand, in pp collisions at GSI, associated 
with the region 0.2 ^ 0.7, the relevant formulae ffTOl) . f|T^ for the asymmetries are 

^ The value of bo/bsp in principle depends on the input values for Q and gNP, but in practice bo/bsp — 1 
GeV, irrespective of Q and gjsip (see [SKTOl for the detail). This in turn implies that (|10p . and thus £^tt{Qt) 
as well as £/'P}^^{Qt), are almost independent of the value of qnp- We have explicitly checked this point 
also in the direct numerical evaluation of ^ with ([2]), ^ in the range qnp = 0.3-0.8 GeV^. 
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Figure 4: The NLL+LO asymmetry £/tt{Qt) of ([7]) with GSI kinematics, 5" = 210 GeV^ 
and y = (f) = 0, and with g^p = O.SGeV^, using the NLO transversity distributions corre- 
sponding to (a) the Soffer bound and (b) the relation (fTTI) . respectively. The corresponding 
values of the saddle-point formula (fTOl) are also shown by the symbols at Qt = 1 GeV. 



dominated by the contributions from the valence components, such that 

^tt{Qt) ^ -c^rr [Qt) ^ - cos(20j „ ^ ^ „ ,2 /;,2 \ ' l^^j 



and ~ [cos{2(f))/2][6u{xl,Q'^)6u{xl,Q'^)/u{xl,Q^)u{xl,Q'^)], using ([12]) and the sim- 
ilar relation for the unpolarized distributions (see (jlj)). As demonstrated in Fig. 3, the 
scale dependence of the w-quark distributions cancels between the numerator and de- 
nominator of (fT6l) in the relevant "valence region". This fact explains why £^tt{Qt) — 
■^tt^{Qt) — -^TT -^iS- 2(b). We note that actually the same mechanism arises for 
all kinematics at GSI, and also for the input transversity distributions using f|TT]) (see 
Fig. 3). This implies that the property, £^tt{Qt) - £^tt^{,Qt) 

~ £/r^^, is characteris- 
tic of all pp collisions at GSI. Moreover, a similar logic applied to ([2]) allows us to derive 
the second approximate equality of (fT5|) : using ( |T2l) and the property 5m(x? 3, ^0/^^) — 
u{xl„ bl/b')6u{xl„ bl/blp)/u{xl„ bl/blp) implied by Fig. 3, (ED gives ArX^LL^ga ^ q2^ ^) 
^ [5u{xl bl/blp)5u{xl bl/blp) /u{xl bl/blp)u{xl bl/blp)]X''^\Ql, Q\ y), up to the 0(a,) 
corrections associated with the coefficient functions (A^) Cf^{z) in ([2]) and ([3]). This imme- 
diately gives the second relation in (fT5!) . combined with ([8]) and (ITOl) . This derivation using 
the properties in the valence region also explains why s^tt{Qt), '^tt^{Qt) in VP collisions 
at GSI are flatter than in pp collisions as noted below (fT5l) (see also Figs. 4-7 below). 

Using the same nonperturbative and kinematical inputs as in Fig. 2, Fig. 4(a) shows the 
dependence of the NLL+LO asymmetry s^tt{Qt) of ([7]) on the dilepton mass Q] the dashed 
curve in Fig. 4(a) is the same as the solid curve in Fig. 2(b), and the results in Fig. 4(a) 
show the "maximally possible" asymmetry, i.e., should be considered as optimistic estimate 
(see the solid curve in Fig. 1). We expect that a more realistic estimate is provided by the 
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S offer bound 


Eq. ( 


mD 


Q 


2.5GeV 


4GeV 


6GeV 


8GeV 


2.5GeV 


4GeV 


6GeV 


8GeV 


SP 
NB 


21.0% 
19.6% 


26.0% 
25.3% 


31.0% 
30.8% 


34.8% 
34.7% 


7.6% 
7.2% 


12.6% 
12.4% 


18.5% 
18.5% 


23.4% 
23.4% 



Table 1: The values of ^^^^{Qt = 0) of ([8]) for the cases of Figs. 1 (a) and (b), as "Soffer 
bound" and "Eq. (|TT1) ". respectively. SP is obtained using the saddle-point formula (fTOl) 
and NB is obtained using the numerical 6-integration of ([2]) and (IHl). 



— Q = 2.5GeV 
--- Q = 4GeV 

- - Q = 6GeV 



"0 0.5 1 1.5 2 

Qt (GeV) 

Figure 5: Same as Fig. 4(b), but for = 80 GeV^. 

results in Fig. 4(b), which is same as Fig. 4(a) but for the case with the input transversity 
distributions using fITT]) corresponding to the dashed curve in Fig. 1. As Q increases, 
■^tt{Qt) increase, preserving the characteristic flat behavior as functions of Qt- The 
results in Fig. 4(b) using (fTT]) are smaller compared with the corresponding Soffer bound 
results in Fig. 4(a), but still yield rather large asymmetries. 

The saddle-point formula (fTOl) has a particularly simple structure, including only the 
parton distribution functions at the fixed values x = and at a single scale bo/bsp — 
1 GeV, and, combined with f|T5l) . allows us to obtain £^tt{Qt) at Qt ~ 1 GeV. In Table 1, 
the row labeled "SP" lists the result^ obtained using the saddle-point formula ffTOl) . with 
the kinematics and nonperturbative inputs of Figs. 4 (a) and (b); for convenience, the 
corresponding values are plotted by the symbol "□" at Qt = 1 GeV in Figs. 4 (a) and (b). 
"NB" in Table 1 hsts £^^^^{Qt = 0) obtained from ([H]) using the numerical 6-integration of 
([2]) and ([3]). We observe that the simple formula (ITOl) has indeed the remarkable accuracy, 
reproducing the results of NB as well as the NLL-I-LO £^tt{Qt) to 10% accuracy, i.e., to the 
canonical size of 0{as) corrections associated with the NLL accuracy (see the discussion 
below ([9])). We also emphasize that the saddle-point formula (ITOl) indicates that the Q 

II These results correspond to "SP-II" discussed in p,[TU]. 
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Figure 6: Same as Fig. 4(b), but for (a) S = 30 GeV^ and (b) S = 45 GeV^. 

dependence of s^tt{,Qt) is controlled by the detailed 2 dependence of the transversity 
and density distributions at the scale ho/hsp. Namely the Q dependence in Fig. 4 can be 
understood by the behavior of the solid lines in Fig. 3 using f[T^ . 

In Fig. 5, we show the NLL+LO asymmetry s^tt{Qt) of O at an another kinematics in 
the collider mode, S" = 80 GeV^ and y = 0, with the transversity distributions using (ITT]) . 
The results are displayed similarly as Fig. 4. The observed pattern is similar as Fig. 4(b), 
but the asymmetries are considerably larger. 

Figure 6 shows the NLL+LO asymmetries s^tt{Qt) in the fixed-target mode, which is 
associated with more challenging kinematic regime for the application of QCD factorization 
framework: (a) and (b) plot the results for 5* = 30 and 45 GeV^, respectively, with y = 
using f|TT]) . The dot-dashed curve in (a) ends at the kinematical upper bound for the 
partonic subprocess, Qr.max = — (a;5)^][l — (a^2)^]/(^i + ^2) — 1-3 GeV; because of 

this kinematical bound, we do not show the cases with Q higher than 4 GeV. Compared 
to (a), the case (b) probes smaller x'j'g region, resulting in the smaller asymmetries. Both 
results give larger asymmetries than the corresponding collider results with the same Q. The 
characteristic behaviors as functions of Qt and Q emphasized above are again observed, 
and actually the specific contributions associated with the results in Fig. 6 obey similar 
pattern as in Figs. 2 (a) and (b), with e.g., £^tt{Qt) — •^TT^iQT) — s^tt- 

We mention the rapidity dependence of the NLL+LO asymmetry ([7]). Figures 7 (a) and 
(b) are same as Figs. 4(b) and 6(b) in the collider and fixed-target modes, respectively, but 
for y = 0.5 and 0.3, and all results agree nicely with the saddle-point results using f fTOj) . 
The characteristic behaviors common to the preceding results are again observed, but the 
values of asymmetry are slightly smaller than the corresponding values for y = 0. 

Figure 8 summarizes our results for the double transverse-spin asymmetry £^tt{Qt) 
of ([7]) at GSI as functions of the dilepton mass Q with y = and various values of 5" 
corresponding to the fixed-target as well as collider mode. The symbols "a" and "v" 
plot s^tt{Qt = QT,pcak) using the NLO transversity distributions corresponding to the 
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Figure 7: (a) Same as Fig. 4(b), but for y = 0.5. (b) Same as Fig. 6(b), but for y = 0.3. 



Soffer bound and ffTTjl . respectively, where QT,peak denotes the value of Qt at which the 
peak of the corresponding NLL+LO tDY cross section is located (see Fig. 2(a)); e.g., the 
several symbols "v" in Figs. 8 (a), (b), (c), and (d) denote the values of the curves at 
Qt = Qr.pcak (— 1 GeV) in Figs. 6(a), 6(b), 5 and 4(b), respectively. The dashed and 
dot-dashed curves draw the results of the saddle-point formula (fTOl) using the same parton 
distributions as those for "a" and "v", respectively. The results remind the readers of 
the relation (fT5l) . Combined with the result in Fig. 3, this relation convinces us that the 
property, ^tt{Qt) — ■^tt^{Qt) — ^tt^ which was discussed below ( IT6l) . indeed holds at 
all relevant kinematics at GSI. 

The conventional double transverse-spin asymmetry Att is defined as the ratio of the 
Qr-integrated polarized and unpolarized DY cross sections, and Att at GSI kinematics has 
been studied in previous works [71 [11], [12]. It is worth noting the relation of Att with our 
Q-r-dependent asymmetry £^tt{.Qt)- As explicitly demonstrated in [10], the integral of ([1]) 
over Qt reproduces exactly the corresponding NLO cross sections for the Q-p-unobserved 
case; indeed A = 0at6 = 0in([H]) ensures [131 [HI E] that the contributions of 0{a1) or 
higher from the resummed components ([2|), ([3|) vanish in the Qt integral. As a result. 



Att = 



/ dQ% {^Td(j/dQ'^dQldyd(t)) 
JdQl {da/dQMQldyd(j)) 



■ cos(2</)) 



^if(x°,x°;g^) + 
ff(x?,x0;g2) + . 



(17) 



where the ellipses stand for the NLO correction terms, coincides completely with the NLO 
asymmetry calculated in [H]. It has been found [H] that the effects of the ellipses in ([T7I) 
on Att are generally small at GSI; i.e., the corresponding K factors of the polarized and 
unpolarized cross sections are similar to each other and cancel out in the ratio. Indeed the 
LO contribution shown explicitly in ([T71) coincides with ([14]), i.e., = while the 

NLO calculation yields e.g., Att = 25.0% for the case of Fig. 2(b). 

In view of the fact Att — ^tt GSI, it is also interesting to calculate the conventional 
LO asymmetry A^j, using the LO transversity distributions, obtained by Anselmino et al. [1] 
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Figure 8: The double transverse-spin asymmetries as functions of Q at y = 0. The dashed 
and dot-dashed curves plot the results of the saddle-point formula ffTOl) . using the NLO 
transversities corresponding to the Soffer bound and flTT]) . respectively, which yield the 
values of ^tt{Qt) of ([7]) for Qt = Qr.pcak (— 1 GeV) as indicated by the triangle up 
and down symbols. The two-dot-dashed curve shows the conventional Qr-independent LO 
asymmetry using the LO transversities corresponding to the Soffer bound. The dotted 
curve shows the "upper bound" of implied by the result of the global fit in [1]. 



through the global fit to the data. (For the unpolarized distributions in the denominator 
we use LO GRV98 distributions [H].) The upper limit of the one-sigma error 
bounds of their fitted transversities (shaded area in Fig. 1) yields the "upper bound" of the 
corresponding as shown by the dotted curves in Fig. 8. Also plotted by the two-dot- 
dashed curves are A^^ using the LO transversity distributions corresponding to the Soffer 
bound, which are given by the dot-dashed curves in Fig. 1. We here note that the above 
mentioned relations, s^tt{Qt) — £^tt^{Qt) — -^tt ^"^^ ^tt = s^tt^ imply ^tt{Qt) — 
Al^. Therefore, the dotted curves in Fig. 8 represent not only estimate of Att-, but also 
that of £^tt{Qt ~ Qr.peak), using the empirical information of the transversities available 
at present, up to anticipated modification of the empirical transversities at the NLO level. 
The effect of the corresponding NLO-level modification of the input transversities on the 
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asymmetries might not be so large, as suggested by the fact that the dashed curve in each 
panel of Fig. 8 is close to the corresponding two-dot-dashed curve. In the small Q region, 
our full result of ^tt{Qt) using f|TT]) can be consistent with estimate using the empirical 
LO transversities, but these results have rather different behavior for increasing Q, which 
reflects the different x-dependence of the corresponding transversities shown in Fig. 1. This 
demonstrates that the experimental data to be observed at GSI, in particular the behavior 
of the Qr-dependent as well as Qr-independent asymmetries as functions of Q, will allow 
us to determine the detailed shape of the transversity distributions in the valence region. 

To summarize, the double transverse-spin asymmetry at a measured Qt to be observed 
at GSI is very useful to determine the transversity, and is complementary to the conventional 
asymmetry associated with the Q^-integrated cross sections: both asymmetries are large 
at GSI, and actually the values of these two asymmetries are almost the same for all GSI 
kinematics. This (approximate) "equality" of the two asymmetries is unexpected from the 
outset and characteristic of pp collisions in GSI experiments, and we have revealed nontrivial 
roles played by the soft-gluon-resummation contributions. Indeed, only after performing 
the soft gluon resummation, we obtain the "physical" behavior for the Qt spectra of the 
DY lepton pair and the associated asymmetry s^tt{Qt) in the small Qt region, with a well- 
developed peak for the former and the flat behavior for the latter, and thus we are able to 
make the reliable estimate. Furthermore, the resummation modifies the parton distributions 
involved in £^tt{Qt) into those with the "effective" scale around Qt ~ 1 GeV, instead of 
Q in the conventional asymmetry. Another reason for the above "equality" is the similarity 
of the QCD evolution between the transversity and unpolarized quark distributions in the 
valence region relevant to the GSI kinematics. These mechanisms have been explicitly 
embodied by the novel saddle-point formula, which relates the asymmetry ^tt{Qt) with 
the transversity distributions at the scale around 1 GeV, in a way as simple as in the 
conventional LO asymmetry. Thus GSI measurements of the asymmetries at a small Qt 
for a variety of dilepton mass Q directly probe the shape of the transversity distributions. 
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